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RECOMB  IN  VTION  AND  IONIZATION  IN  A  NITROGEN  PLASMA 


I.  INTRODUCTION 

When  air  is  heated  and  is  highly  ionized,  nitrogen  and  its 
ions  become  relatively  abundant  for  electron  temperatures  of  1-3 
e.v.  Emissions  from  such  air  provide  diagnostic  information  and 
play  a  role  in  coding  the  heated  channel.  To  describe  such  a 
channel  one  must  therefore  perform  a  radiative  transfer 
calculation  which  requires  a  detailed  knowledge  of  the  air 
species  including  excited  states.  To  provide  this  information, 
we  develop  as  a  first  step,  a  collisional-radiative  model  for 
nitrogen  and  its  ion. 

It  is  well  known  that  a  recombining  plasma  is  a  source  of 
radiation.  The  character  of  this  radiation,  line  versus 
continuum  emission  and  absorption  or  uv  versus  visible  and  infra 
red  radiation,  depends  on  detailed  aspects  of  the  plasma.  In 
addition  to  species  type  important  properties  include  electron 
density,  electron  temperature,  and  optical  character  (thin, 
thick,  or  optically  thick  in  select  spectral  regions).  To 
quantitatively  describe  the  spectral  character  of  radiation 
emitted  from  a  plasma  these  properties  must  be  incorporated  in  a 
model  which  accurately  monitors  the  evolution  of  the  plasma, 
determines  the  extent  to  which  the  plasma  is  or  is  not  in 
equilibrium,  and  accounts  for  the  radiative  transfer  accordingly 
For  example,  in  the  simplest  case  of  a  plasma  in  local 
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thermodynamic  equilibrium  (LTE)  the  radiative  transfer  equation 
may  be  solved  assuming  Planck's  function  correctly  represents  the 
source  radiation  [1].  A  more  complete  description,  however, 
requires  coupling  the  radiative  transfer  equations  to  the  rate 
equations  which  provide  the  population  densities  of  all  bound  and 
continuum  states  involved  in  the  radiative  and  collisions! 
processes  of  the  specific  plasma.  The  work  presented  in  this 
paper  is  the  initial  step  in  constructing  such  a  model. 

A  colli3ional-radiative  model  developed  by  Bates  et .  al.  [2] 
has  provided  the  framework  for  studies  which  give  the  excited 
state  populations  and  their  deviation  from  LTE  for  hydrogen  and 
hydrogenic  ions.  Their  theory  is  statistical  in  nature  and 
treats  the  problem  in  a  simple  and  tractable  way.  Effective 
collisional-radiati ve  recombination  and  ionization  coefficients 
are  defined  and  determined  from  a  knowledge  of  the  rate 
coefficients  for  collisional  excitation  and  ionization  and  the 
spontaneous  transition  probabilities,  for  the  case  of  an 
optically  thin  plasma.  The  influence  of  self-absorption  is  then 
demonstrated  for  the  special  cases  of  a  hydrogen  plasma 
optically  thick  for  various  optical  transitions.  Their  approach 
is  to  solve  a  coupled  set  of  rate  equations  for  the  steady  state 
values  of  excited  state  densities  at  fixed  electron  densities  and 
temperatures,  thereby  showing  the  extent  to  which  the  excited 
state  populations  deviate  from  Saha  equilibrium. 
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The  work  of  3ates  et .  al.  has  motivated  others  to  study  the 
effects  of  coliisional  and  radiative  processes  in  plasmas. 

Orawir.  13 J  has  used  more  recent  cross  section  information  and 
recalculated  the  collisional-radiative  recombination  and 
ionization  coefficients  for  hydrogen  and  hydrogenic  plasmas 
analogous  to  Bates  et.  al . ,  but  within  the  context  of  a  more 
general  formalism.  Park  [4]  has  used  the  method  of  Bates  to 
compute  spectral  line  intensities  for  specific  non-equilibrium 
conditions  in  a  collision-dominant  nitrogen  plasma.  Kulander  [5] 
has  also  studied  the  spectral  properties  of  radiation  emitted 
from  a  non-equilibrium  nitrogen  plasma  as  well  as  the  deviation 
from  Saha  equilibrium  of  the  excited  and  ionized  states  of  an 
optically  thin  plasma. 

In  the  present  work  the  focus  is  on  an  optically  thin 
nitrogen  plasma,  however,  unlike  Bates  and  others  we  solve  the 
time-dependent  coupled  rate  equations  explicitly,  in  a  manner 
similar  to  the  approach  taken  by  Ali  &  Jones  [6]  to  study 
recombination  lasers  in  hydrogen  and  hydrogen  like  plasmas.  The 
calculations  are  for  a  given  electron  temperature,  but  allow  the 
electron  density  to  evolve  in  time.  As  expected  the  population 
densities  show  an  early  transient  behavior  followed  by  a  steady 
state,  from  which  we  calculate  effective  collisional-radiative 
recombination  and  ionization  coefficients  and  compare  the  final 
state  populations  to  those  expected  if  Saha  equilibrium 
conditions  were  satisfied.  Our  rate  equations  include  the  most 


satisfied.  Our  rate  equations  include  the  most  recent 
ionization,  excitation,  and  recombination  rate  coefficients 
available,  some  of  which  are  similiar  to  those  used  by  Drawin  [3 
and  Bates  12],  and  some  are  based  on  experimental  data.  Finally 
we  consider  the  case  of  a  nitrogen  plasma  optically  thick  in 
specific  uv  bound-bound  sDectral  lines. 


The  coupled  rate  equations  which  describe  the  evolution  of 


the  nitrogen  plasma  include  the  following  collisional  and 
radiative  processes: 


Collisional  excitation  -  excitation  of  an  upper 
electronic  state  by  collisions  of  electrons  with  atoms 
(ions)  in  a  lower  state. 

e-N(i)  - >  M ( j )  +  e  j  >  i  X(i,j) 

*  +  +  ( 1 
e*>i(p)->N(q)+a  q  >  p  X  ( p  ,  q ) 


fc; 

t- 


Collisional  de-excitation 
9  *  M  (  j )  ->  M ( i )  +  e 

e  +  N*(q)  ->  N  +  ( p  )  +  e 


-  the  inverse  of  1 . 
j  >  i  Y(j,i) 

q  >  ?  Y+(q,p) 


(2 


3. 


Collisional  ionization  -  ionization 
e  +  N(i)  ->  e  +  e  +  N* ( p ) 
e  +  M+(p)  ->  e  ♦  e  +  N++(u) 


of  atoms  and  ions 
S( i ,p)  ^ 

S+(p,u) 


4.  Three  body  recombination  -  the  inverse  of  3* 
e  e  +  N+(p)  ->  N(i)  +  e  a^(p,i) 

e  +  e  +  >J+*(u)->N*(p)+e  a^Cu.p) 


(4 


5.  Radiative  recombination, 
e  +  N  +  ( p )  ->  N(i)  +  hv 

N+*(u)  ->  N+(p)  +  h v 


aR(p .  i  ) 
oj(u.p) 
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e  + 


(5) 


Spontaneous 

emission. 

N(  J) 

1 

V 

H- 

)  +  h\j 

j  >  i 

A(  j.i) 

N+(q) 

->  N  + ' 

p )  ♦  hv 

q  >  p 

A*  (  q  ,  ? 

Since  we  are  initially  assuming  a  nitrogen  plasma  that  is 
optically  thin  for  all  spectral  regions  photo-ionization  and 
absorption  have  been  excluded.  These  processes  are  accounted  fo 
most  easily  by  describing  photo-ionization  as  negative  radiative 
recombination  and  absorption  as  negative  spontaneous  emission  as 
was  done  by  Drawin  [3],  and  will  be  included  in  detail  in  a 


subsequent  report. 


VO 


rate  equations  and  rate  coefficients 


A.  Coupled  rate  equations 


Let  N"  1  (31)  denote  the  population  density  of  the 
(s-l)-th  ion  (or  neutral  if  2-1)  in  the  electronic  state  whose 
index  is  31,  where  it  is  understood  that  the  density  depends  on 
time.  The  rate  equation  describing  the  evolution  of  this  state 
subject  to  the  coliisional  and  radiative  processes  outlined  in 
the  previous  section  is  given  by 


dN  2_1(m)  =  Me  [  I  N2_1(i)  X2_1U,m)  -  N2-1  (a)  7  Y2"1  (m,  l)  ] 
It  i<m-1  Km  - 1 

+  Me  [  l  N2-1 (n)  Yz"1(n,m)  -  N2'1 (m)  l  X2-1 (m,n)  ] 
n_>m-1  n>m-l 

*  Me  [Ne  7  N2(n)a2_1 (n,m)  ♦  l  N2 ( n ) a2-1 ( n , m) ] 
n  3  n  R 


Ne  N2  1 (m)  l  S2'1 (m,n) 


-  Me  [Me  Nz  1  (m)  7  a 2  2(m,Z)  +•  M2  1  (m)  7  2(m,2,)j 

l  J  i 


Me  I  Mz~2(l)  S2"2  U,m) 


*  l  M 2-1  (n)  A 2 "  1  ( n  ,  m )  -  Nz"'(m)  7  n2"'U)  A2_1a,m: 
n  >  m  1  l  <  m  - 1 


where  We  is  the  electron  density.  The  time  history  of  ail 
species  is  obtained  by  salving  the  set  of  equations  constructed 
from  equation  (7)  for  all  ions  and  neutrals  in  all  states,  whose 
indices  are  n,  L,  and  m. 

3.  Energy  Level  Model 

We  have  chosen  to  study  the  nitrogen  plasma  for 
electron  temperatures  in  the  range  of  1.0  e.v.  to  3.3  e.v..  In 
this  temperature  range  the  important  radiation  is  produced  by 
constituents  which  are  atomic  or  ionic.  In  our  calculation  we 
have  included  the  lowest  thirteen  (13)  levels  for  M  I,  the  lowest 
seventeen  levels  for  M  II,  and  two  effective  levels  for  N  III. 
Details  are  presented  in  Tables  I  through  III;  this  information 
can  be  found  in  Weise  [7],  We  note  that  this  model  includes  as  a 
subset  the  levels  for  M  I  and  N  II  used  in  Xulander's  studies 
15],  as  well  as  most  of  the  N  III  levels.  While  not  as 
comprehensive  in  describing  N  I  as  Park  [4],  we  account  for  M  II 
and  N  III  which  he  fails  to  do,  since  his  interest  is  in  a  region 
of  lower  electron  densities  and  higher  temperatures. 

Within  this  model  the  important  bound -bound  radiative 
transitions  are  listed  in  Table  IV;  each  transition  is  described 
by  its  wavelength,  the  number  given  by  Weise  ]7],  spectral 
character,  transition  index  (  for  example,  7 «■  -  1  denotes 
emission  from  the  seventh  to  the  first  level  in  the  icr.  N  II  ), 
and  transition  energy.  This  table  shows,  for  instance,  that  the 
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visible  radiation  comes  largely  from  transitions  in  tne  ion  NT  I 
Therefore,  one  would  expect  this  radiation  to  be  important  once 
electron  temperatures  become  sufficiently  high  that  these  upper 
levels  of  N  II  are  accessed. 

The  important  free-bound  radiative  transitions  are  listed  in 
Table  V.  These  include  all  transitions  for  the  corresponding 
ions  considered  by  Kulander  [5b J.  While  there  is  a  large  amount 
of  uv  radiation  the  radiative  recombination  from  N+  ->  N(i)  where 
i  _>  7  is  a  source  of  visible  emission.  This  becomes  the 
predominant  source  of  visible  radiation  at  lower  electron 
densities  and  temperatures;  such  that  the  abundance  of  N  II  that 
is  created  resides  in  the  ground  state. 

C.  Hate  Coefficients 

The  rate  coefficients  for  the  collision  and  radiation 
processes  outlined  in  section  II.  and  used  in  the  rate  equations, 
i.e.  equations  (7)  were  obtained  from  several  sources. 

1.  Transitions  amongst  the  levels  of  the  ground  state 

4  2  2 

configurations  of  nitrogen  (  S,  D,  P)  and 
•3  1  ] 

nitrogen  ion  (  ?,  D,  S)  are  optically  forbidden. 
Electron  impact  excitation  coefficients  for  the 
transitions  in  N  I  were  obtained  from  reference 
[3].  Transitions  among  the  metastable  states  in  M 
II  were  accounted  for  using  the  analytic 
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expressions  :or  one  ce-excitation  rate  coefficients 
given  in  reference  [3].  The  reverse  rate 
coefficients  were  obtained  from  detailed  balance. 


2.  Collisional  excitation  rates  of  transitions  which 
are  optically  allowed,  see  Table  IV,  are  obtained 
using  the  following  simple  formula  [9]: 


X ( m , n )  =«  4.3  x  IQ  ^f(n,m)  exp  {-E(n,m)/Te}  ( n , m ) 

E(n,m)  Te1 '  2 


where  S(n,m)  ■  E(n)  -S(m),  m  <  n,  f(n,ra)  is  the 
oscillator  strength,  and 


V  (n,m)  »  0.2  for  ions 


C n , m 5  - ( 1 . 0 * 5 ( n , m ) ) ~1  { (20 .Q+E(n,a) ) >ln[1 .25(1 .0+Te/E(n,m) ) ]} 


for  neutrals 


The  de-excitation  rates  are  calculated  from 


detailed  balance. 


Collisional  ionization  from  a  specific  level  m  to 
the  appropriate  level,  n,  of  the  next  highest  ion 
was  obtained  using  the  formula  of  Drawin  [3], 


a 

■W* 


* 


% 


V 


where  E._, 

is  1 3-59  e . v . , 

5  (m)  the 

number  of 

equivalen 

t  electrons,  AE 

-  Ez(n;  - 

Ez  '■  m )  ,  f  ( m 

the  oscil 

later  strength, 

and  z  *  1 

for  neutrals 

>  1  for  ions.  "  is  given  by 


'•  ‘  AS  i 

Te  ’ 


dx  |1-AE;  ex?  ( - x )  in  ■  .253x 
at/Te  Te  j  AE/Te  j 


(11) 


where  3  =  1 +[ ( z -1 ) / ( z +2  )  ] .  The  oscillator  strength 
was  chosen  so  that  rate  coefficients  given  by 
equation  (11)  matched  closely  the  ionization  data 
in  Table  VIII  of  reference  1 10]  which  is  based  on 
experimental  cross  sections  for  the  ionization  of 
nitrogen;  our  choice  was  f(m)  »  0.2.  The  rate 
coefficient  for  three-body  recombination,  obtained 
from  detailed  balance,  is 


z-1  ,  , 

(n,m) 


1.66  *  10~22  gz  '  (m)  Sz  '  (m,n)  exp  \  AS} 


.  - 1 


z-1 


gZ(n) 


Te 


3/2 


l  TeJ 


(12) 


2  ■”  1 

where  g  (m)  is  the  statistical  weight  of  state  m 
and  gz(n)  is  the  weight  of  the  ionic  state. 

4.  The  rate  coefficients  for  radiative  recombination 
were  obtained  using  the  approximate  formula  of 
Bates  [ 2 ] ,  e . g . 


a*'1!?,!)  ■  5.2 


In-14  .  r._>3/2 

*  1  0  exp  |  Ac.  1  i  A_E 

l  Te  j  [ Te  J 


E  f  AE 
:  Te 


Here  E,  is  the  exponential  integral 


s,  fiEl 

(Tel 


dx  exp  fAE  x 
x  i  Te 


(13) 


(14) 


.v  .•  . ••  .  .  - 


-Vu.  -1 . 


•  v  '  -  >‘v1 

.M  .  -  N 

Oa.aU.'.v'-V.vj 
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F or  example,  eq.'.  ’3)  gives  the  rate  coefficient  for 
reccmb i nat i on  of  level  p  of  N  II  to  level  i  of  M  I, 
see  II.  5  .  ,  if  z  -  1. 

5 .  Finally,  the  Einstein  A  coefficients  for 

spontaneous  emission  were  taken  from  Weise  [7]. 

3.  Ionization  Potential  Reduction 

The  ionization  energy  of  an  isolated  atom  or  ion  is 
different  from  one  immersed  in  a  plasma  [11]  and  must  be 
corrected  accordingly.  The  net  ionization  potential  reduction  is 
due  to  energy  released  upon  placing  an  electron-ion  pair  in  a 
plasma;  one  of  the  consequences  is  that  electrons  occupying 
states  which  in  an  isolated  atom  would  be  bound  now  become 
unbound  should  the  corrected  ionization  potential  be  less  than 
the  energy  of  the  state.  Our  rate  equations  have  been  modified 
to  account  for  this  high  electron  density  (low  electron 
temperature)  phenomenon. 

At  each  step  in  the  numerical  integration  of  the  rate 
equations,  eqs.  (I7),  the  electron  density  and  temperature  are 
noted  and  the  ionization  ootential  reduction,  ,  calculated 

*  CD  ' 


*w  _  *»  3 

where  AZ“  ,  end  Te  are  in  e.v.  and  Me  has  units  »A/cm  .  Reducing 

the  ionization  potential  an  amount  determined  by  eq.  (15;  has 

the  effect  of  modifying  the  coliisior.al  ionization  rate 

ooefficient,  eq.  (10),  as  follows: 


Sz"  ' (m,n )  »  3Z~ ' (m,n)  x  exp { AEZ“ ' /Te } 


(16) 


For  example,  in  the  extreme  case  of  low  temperature,  Te  -  1.0 

19  3 

e.v.,  and  high  electron  density,  Me  ■  10  /cm  ,  the  ionization 
rats  increases  by  39^  so  its  effect  can  be  substantial. 
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NUMERICAL  RESULTS 


T  *  f 
4 

A.  Time-Dependent  Population  Densities 

Figures  (1)-  (10)  show  representative  results  obtained 

by  integrating  equations  (7).  Numerical  integration  was 

accomplished  by  using  CHEMEQ,  a  routine  specifically  designed  to 

solve  stiff  ordinary  differential  equations  [12].  The  initial 

conditions  for  this  set  of  figures  are  Te  »  2.0  e.v.,  N(1)  -  4. 

x  10’^/cm^,  Me  =  N  +  (1)  =  1  .0  x  lO^/cm^  (i.e.  charge  neutrality) 

with  all  other  population  densities  being  set  equal  to  zero 

4  R 

(.actually  1  x  10  /cm  )  .  The  general  behavior  of  the  time- 
dependent  solutions  is  as  expected  -  an  early  transient  behavior 
followed  by  rapid  relaxation  to  a  steady  state.  The  relaxation 
time  is  on  the  order  of  15  -  22  nsec.  Such  rapid  relaxation  led 
3ates  et .  al.  [2]  to  solve  for  the  excited  state  population 
densities  using  a  steady  state  approximation. 

In  Figure  (1)  the  total  population  densities  of  the 
electrons,  neutrals,  and  ions  are  presented.  At  Te  *  2.0  e.v. 
the  number  of  N++  ions  is  negligible,  i.e.  approximately  three 
orders  of  magnitude  less  than  the  neutrals.  Consequently,  the 
total  number  of  electrons  nearly  equals  the  total  number  of  N+ 
ions  which,  in  the  steady  state,  exceeds  the  population  of 
neutrals  by  a  factor  of  three. 


r  v  •.  v  v  vv.  ’rvir  twv  "  - » - 

Examination  of  Figures  (2)  -  (4)  shows  that  most  of  the 
nitrogen  population  resides  in  the  ground  state  and  the 
metastable  states,  particularly  the  20°  state.  The  levels  which 
contribute  to  the  emission  of  bound -bound  uv  radiation, 

N(4-6,’3),  are  approximately  equal  in  population  to  those  which 
contribute  to  the  emission  of  bound-bound  ir  radiation,  JJ(9-12). 

As  in  the  case  of  N  the  ground  and  metastable  states  of 
N  +  are  the  more  heavily  populated,  see  Figures  (5)  -  (7).  Unlike 
M  only  M  (4)  and  N+(5),  i.e.  the  first  two  excited  states  which 
emit  bound-bound  uv ,  are  highly  populated  because  of  the  low 
electron  temperature.  The  rest  of  the  uv  and  visible  emitting 
levels  are  two  orders  of  magnitude  less  than  the  N  radiating 
levels.  The  population  of  the  lower  N  +  levels  is  important  for 
the  emmision  of  radiation  resulting  from  radiative  recombination, 
f  ree-bound . 


Figure  (8)  shows  that  the  ions  that  are  created 

reside,  as  expected,  predominately  in  the  ground  state. 

While  a  direct  comparison  cannot  be  made  with  the 

results  of  Kulander  [5a]  our  steady  state  populations  are  in 

qualitative  agreement.  The  relative  population  of  the  major 

1 6  7  1  Q  7 

equilibrium  species  for  Ne  -  10  /cm  to  10  /era  and  Te  *  1.0 
e.v.  to  2.0  e.v.  follows  the  trend  obtained  by  Kulander.  At  Te  - 
1.0  e.v.  N(l)  is  the  major  species  for  Me  =■  lO^/cm^  to  IG^/cm^ 
while  at  Te  -  2.0  e.v.  N+(l)  predominates.  At  Te  *  1.5  e.v. 
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the  major  species  at  lower  electron  densities 


/cm"  - 


:nr 


while  M(l)  is  at  higher 


Kulancer  also  presents  results  on  the  distribu 
amongst  species  versus  electron  concentration, 
are  in  general  agreement.  Detailed  aspects  of 
distribution  are  discussed  in  section  IV.  C.. 


/  ’9  7 , 

density  (’O'  /cm" ) . 
ticn  of  population 
Again  our  results 
our  population 


3.  Coilisional-P.adiative  Recombination-Ionization 
Coefficients 

Bates  et .  al.  [2]  have  defined  effective  recombination 
and  ionization  coefficients  which  account  for  the  net  effect  of 
collisions!  ionization,  three-body  recombination,  and  radiative 
recombination  in  a  very  simple  way.  These  coefficients  have  been 
tabulated  exclusively  for  hydrogen  and  hydrogenic  ion  plasmas, 

[2]  and  C 3 ] •  In  this  section  we  define  the  effective 
coefficients  analogous  to  Bates  et.  al.  and  present  results  for 
our  optically  thin  nitrogen  plasma  as  well  as  the  case  where  all 
bound-bound  uv  lines  are  assumed  to  be  optically  thick. 

Denoting  a„a,  and  as  the  effective  recombination 

coefficients  for  M+->  M  and  M**-  ->  M  +  and  and  S*R  as  the 

effective  ionization  coefficients  for  the  reverse  processes  the 
rate  equations  for  N,  N+,  N++,  and  Me  may  be  written 


dM  ’ 


Me  N*  a° 


Me  M  S° 


+  Me  M 


dM*  -  -Me  N* 

—  CR 

—  u 

-  Me 

M  3°R  *  Me 

+  4*  +  4- 

dM  -  Me  M  SCR 

-  Me 

+  4*  4- 

N  3CR 

3 Me  -  Me  M  3°c  ♦ 

i  “  *“  Wlk 

•1  O 

Me  M 

+  S^R  -  Me  M 

Therefore,  in  terms  o 

f  the 

detailed  ra 

section  -I.  a„R,  aCR, 

S° 

CR  ’ 

and  S*R  are 

-  Me  N  +  S„D 

u  a  on 


(17) 


+  4  '  Me  M"  +  aCR 

te  coefficients  outlined  in 
given  by 


„o 

^CR 


{Me  l  a,(p,i)  N  +  (p)  l  a_(p,i)  N+(p)}/M  + 

n  4  j  n  ’• 

9  >  -  P  »  «*■ 

{Me  i  aR(u,p)  N**(u)  ♦  l  a *  (u,p)  M++(u)}/M+* 
p  ,u  J  p  ,  i 

{  l  M(i)  S(i,p)}/N 
i  >  P 

{  l  M*(p)  S+(p,u)}/N+ 

P  »u 


(18) 


A  quick  check  shows  that  these  definitions  satisfy  the  conditions 
for  particle  and  charge  conservation. 


The  coefficients  defined  in  equations  (17)  and  (18)  are 
time-dependent,  due  to  the  time  dependence  of  the  contributing 
population  densities;  see  Figs.  (9)  -  (10).  The  steady  state 
values  are  presented  in  Tables  '/I  -  IX  as  a  function  of  the  final 
electron  density  and  electron  temperature  (Te  was  fixed 
throughout  the  integration).  Close  examination  of  these  results 
reveals  the  following:  At  high  electron  densities  the  effective 
recombination  coefficient  is  linearly  proportional  to  Me.  This 
is  due  to  the  dominance  of  three-body  recombination  over 

1  3 

radiative  (two-body)  recombination  in  this  regime,  Me  _>  1 0 
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1 9  < 

10  /cm".  For  lower  electron  densities  the  opposite  is  true  and 

the  effective  coefficient  approaches  the  two-body  value.  The 

collisional-radiative  ionization  coefficients  show  a  rapid  rise 

1 5  1 7  9 

for  low  electron  densities  (Me  <_  1 0  -  10  '  /cm  )  and  a  much 

1 8  9 

slower  rise  for  higer  densities  (Me  _>  10  /cm  )  .  This  is  in 

contrast  to  the  behavior  of  the  coefficients  of  Dr aw in  and  Bates 

which  show  a  saturation  at  high  Me.  The  density  dependence  in 
o  + 

and  at  high  Me  in  our  calculation  comes  from  the 

correction  discussed  in  section  III.  D..  The  level-specific 
coilisional  ionization  coefficients  depend  on  Me  and  Te  in  the 


manner  of  eq.  ( 1 6 )  ;  this  carries  through  to  the  definition  of  S 
and  in  eq.  (13).  Neither  Drawin  nor  Bates  included  this 

correction . 


CR 


The  collisional-radiative  coefficients  in  Tables  VI  and 

VIII  can  be  loosely  compared  to  the  results  of  Bates  .et.  al .  [2] 

and  Drawin  [3j.  While  there  i3  no  apriori  expectation  that  the 

results  should  be  in  strict  agreement,  to  the  extent  that  a 

nitrogen  plasma  differs  from  a  hydrogenic  plasma,  many  of  the 

rate  coefficients  were  obtained  using  similiar  (hydrogenic) 

expressions.  However,  it  must  be  pointed  out  that  the  results  of 

Bates  et.  al.  and  Drawin  are  not,  themselves,  in  agreement. 

Specifically,  the  recombination  coefficients  of  Bates  et .  al. 

exceed  those  of  Drawin  by  up  to  a  factor  of  four  while  the 

ionization  coefficients  are  greater  by  an  order  of  magnitude.  At 

1 6  3 

low  electron  densities,  Me-  10  /cm,  our  re combi  nation 
coefficients  are  generally  a  factor  of  one  to  four  less  than 
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those  of  Drawin  ant  Bates  et.  al.,  depending  on  temperature.  At 


higher  Me,  -  ’Q^/cm", 


■e combi nat ion  coefficients  fail 


between  their  results.  Cur  col  1  i s ior.al -radi at i ve  ionization 

coefficients  are  in  good  agreement  (factor  of  1.0  -  2.0)  with 

1 6  7  17  7 

those  of  Drawin  at  low  Ne,  10  /cm  to  1 C  /cm".  At  higher  Ne 
our  coefficients  exceed  those  of  Drawin  by  up  to  a  factor  of 
ten.  ?hi3  difference  is  due  primarily  to  the  ionization 
potential  reduction  correction  discussed  above.  As  a  reminder, 
the  calculations  of  Bates  and  Drawin  are  for  hydrogen  where  at 
higher  electron  densities  the  character  of  the  ion  does  not 
change.  In  contrast,  in  nitrogen  other  stages  of  ionization 
arise. 


In  a  subsequent  report  it  will  be  shown  that  the 

absorption  coefficient  for  our  model  nitrogen  plasma  is  on  the 
4 

order  of  10  /cm  in  the  uv  spectral  regime.  This  corresponds  to  a 

-4 

mean  free  path  of  -  10  cm.  In  Tables  X  through  XIII  the 

collisional-radiative  coefficients  are  presented  for  the  case 

where  all  bound-bound  uv  lines  are  assumed  to  be  completely 

reabsorbed.  As  expected,  there  is  very  little  difference  between 

the  optically  thin  and  optically  thick  recombination 

1 3  7 

coefficients.  For  high  electron  densities,  _>  10  /cm  ,  where  the 
plasma  is  collision-dominant  (see  the  discussion  in  the  next 
section),  the  ionization  coefficients  are  comparable.  However, 
for  lower  electron  densities  the  optically  thick  ionization 
coefficients  are  much  larger  than  the  optically  thin  case,  as 
great  as  a  factor  of  ten. 
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When  collisional  transition  nates  exceed 


the 

co r responding  radiative  transitions  the  plasma  is  said  to  be 
collision-dominated.  It  is  then  possible  to  specify  the 
populations  of  the  constituent  species  by  using  the  Boltzmann- 
Saha  equation  at  the  local  electron  temperature  111].  The  onset 
of  Saha  equilibrium  depends  on  Me  and  Te;  for  example,  an 
optically  thin  nitrogen  gas  at  Te  =  1.0  e.v.  becomes  collision- 
dominant  at  Me  -  10l3/cm3  [5a]. 

One  form  of  the  Saha  equation  is,  [11], 


Me  Mz(p) 
MZ~'  (i) 


2  g2(p) 

z-1  ,  .  . 
g  ( i ) 


r 


mkTe 


'STr-n; 


e]  3  72  exp { -[ Ez ( p )  -  Ez ” 1 ( i ) ] /kTe } 


(19) 


We  have  used  eq.  (19),  after  correcting  the  exponential  factor 
for  the  ionization  potential  reduction  (see  section  III.  D.),  to 
compare  the  steady  state  population  densities  with  those  expected 
assuming  a  collision -dominant  environment  for  a  range  of  electron 


dens i t i es , 

Me  = 

[1 

O 

O'* 

-  10 

,9/=m3!. 

and  temperatures, 

Te  =  [1.0 

The 

ra 

t  io 

of  t 

he  level 

popul a t i on  to  the 

Saha 

populat ion 

may 

be 

de  f 

i  ned 

as 

w here  '  J  “ 
s 


( r. '  is  the  solution  to  the  corrected  eq.  (19,'.  In 

Table  XIV  we  present  threshold  values  for  Me  above  which 
_  * 

,n  =  -  ♦_  .-5  for  the  indicated  levels,  i.e.  Saha  equilibriu 
is  achieved  to  within  15*.  Column  1  requires  this  criteria  to  b 
satisfied  for  all  if  I  and  M  II  levels,  column  2  restricts  the 
criteria  to  excited  states  of  N  I  (specifically,  i  -  7  to  12)  > 
and  column  3  bo  excited  states  of  N  II  (i.e.  p  *  7  to  17).  The 
results  show  the  following:  First,  for  any  given  Te,  as  the 
number  of  electrons  increases  equilibrium  is  achieved  for  excite 
M  I  states  before  excited  M  II  states  and,  as  expected,  before 
the  lower  states.  In  addition  to  Me  the  threshold  values  depend 
on  Te.  As  Te  increases  more  electrons  are  neccessary  to  have  a 
completely  collision-dominant  environment,  i.e.  complete  LIE. 

For  an  optically  thin  nitrogen  plasma  a  state  of  complete  TIE 
exists  if  Me  *  (2.2  -  5.0)  x  10’^/cm^  for  Te  ranging  from  1.0  to 
3-0  e.v..  The  criteria  given  by  Griem  [11]  and  others  on  the 
validity  of  LTS  is  that  the  electron  deexcitation  rate  of  the 
resonance  line  must  exceed  ten  times  the  radiative  decay,  i.e. 
for  the  nitrogen  atom  one  must  have 


( d ,  1  )  Me  >  10  A( 4 , 1  ) 


For  Te  =  1.0  e.v.,  Eq .  (21)  gives  Me  ^  2.3  x  10‘®/cm^  in  our 

case.  However,  while  radiative  processes  may  be  important  to  th 
lower  levels  the  excited  states  of  M  I  are  still  controlled  by 
oollisional  processes  for  Me  _>  (5.7  -  2.1)  x  10'  /cm  ever  the 
same  temperature  range.  This  quas i -stat ionary  state  exists  for 


’equ :  re; 


The  recombination  and  ionization  processes  for  an  optically 

thin,  homogeneous,  charge  neutral  nitrogen  plasma  have  been 
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studied  for  electron  densities  between  10  /cm  and  10  '/cm  and 
electron  temceratures  from  1.0  .ev.  to  3-0  e.v..  Effective 


collisional-radiative  recombination  and  ionization  coefficients 


are  calculated  which  describe  the  net  recombination  and 

ionization  and  may  be  used  in  sophisticated  plasma  chemistry 

simulation  codes.  The  bound-bound  and  free-bound  radiation  has 

been  characterized.  It  has  been  shown  that  for  electron 

1 8  3 

densities  greater  than  -  10  /cm  a  state  of  local  thermodynamic 

equilibrium  exists  whereas  exited  states  may  be  in  equilibrium, 

1 7 

collision-dominant,  for  electron  densities  as  low  as  -  10  /cm-'. 

This  has  implications  for  simplifying  the  description  of 

1 3  3 

radiative  transfer  in  a  plasma.  For  example,  above  Me  -  10  /cm 

a  transfer  theory  based  on  the  assumption  of  ITS  should  be 

1 7  3 

adequate  while  below  Me  -  10  /cm  individual  species  must  be 
considered  in  detail.  In  the  intermediate  regime  hybrid  theories 


should 


prove  useful. 
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2.0 

1.3 

1  .02  +  17 

6.6 

4.3 

3.1 

2.2 

1  .5 

2.0 

1  -IE-11 

7.4 

5.6 

3.9 

2.6 

3.0 

1  .7 

1.12-11 

8.1 

5.3 

3.8 

4.0 

2.2 

1  .4 

1  .  0E-1  1 

7.5 

5.0 

5.0 

2.7 

1  .7 

1  -3 

9.4 

6.3 

6.0 

3.2 

2.0 

1.5 

1.12-11 

7.6 

7.0 

3.8 

2.4 

1  .7 

1  .3 

8.9 

3.0 

4.3 

2.6 

2.0 

1  .5 

1  .  0  E  -  1  1 

9.0 

4.8 

3.0 

2.2 

1  .7 

1  .1 

1  .OE+13 

5.4 

3.4 

2.5 

1.9 

1.3 

2.0 

1  -IE-10 

6.8 

4.8 

3.3 

2.7 

3.0 

1  .6 

1 .02-10 

7.2 

5.8 

4.1 

4.0 

2.2 

1  .4 

9.6 

7.7 

5.6 

5.0 

2.7 

1  .7 

1  .22-10 

9.6 

7.0 

6 . 0 

3.3 

2 . 1 

1  .5 

1  .12-10 

8.6 

7.0 

3.3 

2 . 4 

1  .7 

1  .3 

1  .02-10 

3.0 

4.4 

2.3 

2.0 

1  . 5 

1  .2 

9.0 

4.9 

3-2 

2.2 

1  .7 

1.3 

1  .02*1  9 

5.5 

3.5 

2.5 

1  .9 

1  .5 
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<\i  m  rt  in v o  t -  a?  ctn  rvi  ro.tr  in^o  r-co  ^  r-  c\i  rn^t  in  vo  r-  m  o> 


Coilisicnai-Radiative  Recombination  Coefficients 
Optically  Thin  (All  Wavelengths) 

Ion 


e  ( e  .  v  . ) 

Me 

1  .0 

1  .5 

2.0 

2.5 

3-0 

.05+1 5 

3.65-12 

2.1  5-12 

1  .  7E-1  2 

1  .  3  E  - 1  2 

9.75-13 

r*. 

• 

4 . 6 

2.6 

2.0 

1 

1.15-12 

.0 

5.3 

3.1 

2.4 

1  .8 

i  .3 

.0 

5.9 

3.6 

2.8 

2.1 

1  .4 

.  0 

6.4 

4.1 

3.1 

2.3 

1  .6 

.0 

6.9 

4 . 5 

3.5 

2.5 

1  .7 

.0 

7 . 4 

5.0 

3-3 

2.8 

’  .9 

.0 

3 . 0 

5.5 

4.2 

3.0 

2.0 

.0 

8.5 

6.0 

4 . 5 

3.2 

2.2 

.05+17 

9.0 

6.4 

4.8 

3-5 

2.3 

.0 

1  . 5  E- 1 1 

1  .  0 E- 1 1 

3.1 

5.3 

4.0 

.0 

2.2 

1  .5 

1 . 1 E-1 1 

3.1 

5.3 

.0 

2.9 

1  .9 

1  .5 

1.15-11 

7.4 

.0 

3.5 

2.3 

1  .3 

1  .3 

9.2 

.0 

4.2 

2.7 

2.1 

1.6 

1 . 1 E-1 1 

n 

«  V 

4.9 

3.2 

2.4 

1  .9 

1  .3 

.0 

5.5 

3.6 

2.8 

2.1 

1  .5 

.0 

6.2 

4.0 

3.1 

2.4 

1  .7 

.05+1  3 

6.8 

4.5 

3.4 

2.6 

1  .9 

.0 

1 .4E-10 

8.9 

6.4 

5.2 

3.9 

.0 

2.1 

1 . 3E-10 

9.3 

7.9 

5.9 

.0 

2.3 

1  .8 

1  . 3  E  —  1  0 

1  .IE-10 

7.9 

.0 

3.5 

2.2 

1  .7 

1.3 

9.9 

.0 

4.2 

2.7 

2.0 

1  .6 

1  .  2E-1  0 

.0 

4.9 

3.1 

2.3 

1.8 

1  .4 

.0 

5 . 5 

3.6 

2.7 

2.1 

1  .  6 

•  J 

6.2 

4.0 

3.0 

2.3 

1  .3 

.05*1  9 

6.9 

4.5 

3.3 

2.6 

2.0 

l 
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in  vo  t  -  on  a\ 


Coilisional-Raiiative  Recombination  Coefficients 


Optically  Thin  (All  Wavelengths) 
Neutral 


T  e  ( e  .  v  . ) 

Me 

1  .0 

1  .5 

2.0 

2.5 

3.0 

i  .OE  +  1  5 

4 . 1 S-1 4 

4.0E-1 2 

2. 2E-1 1 

2.2E-10 

4.5E-10 

2  -  0 

5.4 

5.8 

7.0 

2.8 

6 . 0 

3.0 

5 . 4 

7.6 

8.3 

3.5 

7.4 

4.0 

7.2 

9.2 

1  .  0E-1  0 

4.1 

3.6 

5 . 0 

7.9 

1  . 1  E-1 1 

1  .2 

4.6 

9.3 

5.0 

8 . 5 

1  .2 

1  .3 

5.1 

1  .IE-09 

7.0 

9.2 

1  .3 

1  .5 

5.6 

1  .2 

3.0 

9.7 

1  .4 

1  .6 

6.0 

1  .3 

9.0 

1  .  0E-1  3 

1  .5 

1  .7 

6.4 

1  .4 

1  .0E+1  7 

1  . 1 

1  .6 

1  .8 

6.7 

1  .4 

2.0 

1  .5 

2.1 

2.4 

9.0 

2.0 

3.0 

1  .9 

2.4 

2.3 

1 . 0E-09 

2.3 

4.0 

2.1 

2.7 

3.0 

1  . 1 

2.5 

5.0 

2.3 

2.9 

3.2 

1  .2 

2.7 

6 . 0 

2 . 5 

3.1 

3.3 

1.3 

2.9 

7.0 

2.6 

3.2 

3.4 

1  .3 

3.0 

3.0 

2.7 

3.3 

3.5 

1.4 

3.1 

9.0 

2.8 

3.4 

3.6 

1  .4 

3.1 

1  .0E+1  3 

2.9 

3.5 

3.7 

1.4 

3.2 

2.0 

3.5 

4.0 

4.1 

1  .5 

3.6 

3-0 

4.0 

4.2 

4.3 

1.6 

3.7 

4.0 

4.3 

4.4 

4.4 

1  .6 

3.8 

5.0 

4.5 

4.5 

4.5 

1  .7 

3.9 

TABLE  IX 

Coiiisional-Radiative  Ionization  Coefficients 
Optically  Thin  (All  Wavelengths) 


Ion 


Te ( e . v . ) 

Ne 

1  .0 

1  .5 

2.0 

2.5 

3.0 

1  . 0  E  + 1 6 

2.9E-21 

4.0E-1 7 

5 . 8  E— 1 5 

1  . 1 E-1 3 

8.4E-1 3 

2.0 

3.9 

5.5 

6.3 

1  .5 

1  .  0E-1  2 

3.0 

4.9 

7.2 

1  .  0  E- 1  4 

1  .3 

1  .2 

4.0 

5.3 

9.2 

1  .3 

2.2 

1  .4 

5.0 

6.6 

1  . 1  E-16 

1  .5 

2.6 

1  .6 

6.0 

7.4 

1  .3 

1  .8 

3.0 

1.8 

7.0 

3.1 

1  .5 

2.0 

3.4 

2.1 

3.0 

8.3 

1  .6 

2.2 

3.7 

2.3 

9.0 

9.5 

1.8 

2.4 

4.1 

2.5 

1  . 0  E+ 1  7 

1  .0E-20 

1  .9 

2.6 

4.4 

2.7 

2.0 

1  .7 

3.0 

4.3 

7.4 

4.6 

3.0 

2.2 

3.9 

5.6 

1  .OE-12 

6.2 

4.0 

2.6 

4.6 

6.6 

1  .2 

7.5 

5.0 

2.9 

5.2 

7.4 

1  .4 

8.6 

6.0 

3.2 

5.7 

8.1 

1  .5 

9.6 

7.0 

3.4 

6.1 

8.6 

1  .6 

1  .  0E- 1  1 

8.0 

3.6 

6.4 

9.2 

1  .7 

1  . 1 

9.0 

3.8 

6.8 

9.6 

1  .3 

1  .2 

1  . 0  E+ 1 3 

4.0 

7.1 

1  . 0  E  —  1  3 

1  .9 

1  .  2 

2.0 

5.4 

9.1 

1  .2 

2.4 

1  .6 

3-0 

6.5 

1  .  0  E  —  1  5 

1  .4 

2.6 

1.8 

4.0 

7.5 

1  . 1 

1  .5 

2.3 

1  .9 

5.0 

3.4 

1  .2 

1  .5 

2.9 

2.0 

6 . 0 

9.0 

1.3 

1  .6 

3.0 

2 . 1 

7.0 

9.5 

1.3 

1  .7 

3.1 

2 . 1 

8.0 

1  .  0E- 1 9 

1  .4 

1  .7 

3.2 

2 ’.2 

9.0 

1  . 1 

1  .4 

1.8 

3.3 

2.2 

1  .02+1  9 

1  .2 

1  . 5 

1  .8 

3.3 

2.3 

in  '-O  c~ 


~o*l i3ionai-Radiative  Recombination  Coefficients 
Optically  Thick  (All  UV  Bound-Bound  Wavelengths) 

Neutral 


fe  ( e  .  v  . ) 
Ne 

1  .  0  E+ 1  6 


1  .02*1  7 
2.0 
3-0 
4.0 
5.0 
6.0 
7.0 
8.0 
9.0 

1  .  0  E  + 1  8 
2.0 
3.0 
4.0 

.0 


1  .OS-1  9 


1  .0 

1  .5 

2.0 

2.5 

3-0 

2.6E-12 

1  .IE-12 

3 . 2E-1  3 

5 . 4  E  - 1  3 

4. 8E-1  3 

2.9 

1  .4 

1  . 1  E-1 2 

6.7 

5.3 

3-3 

1.8 

1  .3 

3.0 

7.0 

3-3 

2.2 

1  .6 

9.4 

8.2 

4.3 

2 . 6 

1.8 

1  . 1  E-1 2 

9.4 

4.7 

2.9 

2.0 

1  .3 

1  . 1  E-1  2 

5.1 

3.3 

2.3 

1  .4 

1  .2 

5 . 6 

3-7 

2.5 

1  .5 

1  .3 

6.1 

4.0 

2.3 

1  .8 

1  .4 

6.6 

4.3 

3.0 

1  .9 

1  .5 

1  . 1 E-1  1 

7.4 

5.6 

3.6 

2.7 

1.7 

1  . 1  E-1 1 

8.1 

5.4 

4.0 

2.2 

1  .4 

1  . 0  E- 1  1 

7.3 

5.2 

2.7 

1  .7 

1  .3 

9.1 

6.5 

3.2 

2.0 

1  . 5 

1  . 1  E-1 1 

7.8 

3.8 

2.4 

1  .7 

1  .3 

9.0 

4.3 

2.6 

2.0 

1  .5 

1 . 0E- 1 1 

4.3 

3.0 

2.2 

1  .7 

1  . 1 

5.4 

3.4 

2.5 

1  .9 

1.3 

1  .IE-10 

6.8 

4.3 

3.3 

2.7 

1  .6 

1  .0E-10 

7.2 

5.8 

4.1 

2.2 

1  .4 

9.6 

7.7 

5.6 

2.7 

1  .7 

1  .  2  E- 1  0 

9.6 

7.0 

3.3 

2.1 

1  .5 

1  .IE-10 

8.6 

3-3 

2.4 

i  .7 

1.3 

1  . 0  E- 1 0 

4.4 

2.3 

2.0 

1  .5 

1  .2 

4.9 

3.2 

2.2 

1  .7 

1.3 

5.5 

3.5 

2.5 

1  .9 

1  .5 

XI 


TABLE 

Collisional-Radiative  Recombination  Coefficients 
Optically  Thick  (All  3ound-3ound  UV  Wavelengths) 


Ion 


Te (e  .  v  . ) 

Ne 

1  .0 

1  .5 

2.0 

2 . 5 

3.0 

1  .  0  E+  1  6 

4.4E-1 2 

2. 1 E-1 2 

1  .  7E- 1 2 

1  .  1 E-1  2 

1  .  0  E- 1  2 

2.0 

4.6 

2.4 

1  .9 

1  .3 

1  .2 

3.0 

5.1 

3-0 

2.4 

1  .5 

1  .4 

4  .0 

5.6 

3.6 

2.3 

1  .7 

1  .6 

5.0 

6.2 

4.1 

3.1 

1.9 

1  .7 

6.0 

6.8 

4.6 

3.5 

2.1 

1  .9 

7 .0 

7.4 

5.0 

3.3 

2.4 

2.1 

3.0 

3.0 

5.5 

4.2 

2.6 

2.3 

9.0 

8 . 5 

6.0 

4.5 

2.3 

2.4 

1  . 0  E+ 1 7 

9.0 

6.4 

4.8 

3.0 

2.6 

2.0 

1 .52-1 1 

1  .  0E-1  1 

3.1 

6.3 

5.7 

3.0 

2.2 

1  .5 

1  .IE-11 

7.7 

6.2 

4.0 

2.9 

1.9 

1.5 

1 .0E-1 1 

8.0 

5.0 

3.5 

2.3 

1  .3 

1  .3 

9.6 

6 . 0 

4.2 

2.7 

2.1 

1.6 

1  . 1  E-1 1 

7.0 

4.9 

3.2 

2.4 

1  .3 

1  .3 

3.0 

5.5 

3.6 

2.3 

2.1 

1  .5 

9.0 

6.2 

4.0 

3.1 

2.4 

1  .7 

1  .  0  E  + 1  3 

6.3 

4 . 5 

3.4 

2.6 

1.9 

2.0 

1  .  4  E  —  1  0 

8.9 

6.4 

5.2 

3.9 

3.0 

2.1 

1 . 3E-10 

9.3 

7.9 

5.9 

4.0 

2.3 

1  .3 

1  .  3  E  - 1  0 

1 .IE-10 

7.9 

5.0 

3.5 

2.2 

1  .7 

1  .3 

9.9 

/*  fs 

O  .  J 

4.2 

2.7 

2.0 

1  . 6 

1  .  2E- 1 0 

■> 

.  j 

4.9 

3.1 

2.3 

1  .  8 

1  .  4 

3.0 

5  . 5 

3.6 

2.7 

2 . 1 

1  .6 

9.0 

6.2 

4.0 

3.0 

2.3 

1  .8 

1  .  0  E  +  1  9 

6.9 

4.5 

3.3 

2.6 

2.0 

o  o  o  o  o  o  o  o  o 


Coll isior.al-Radiati ve  Ionization  Coefficients 


Optically  Thick  (All  Bound-Bound  UV  Wavelengths) 


Neutral 


e  .  v  . ) 

Ne 

1  .0 

1  .5 

2.0 

2.5 

3.0 

E+ 1  5 

3.5E-1 3 

2.0E-1 1 

1  .  6  E  - 1  0 

4. 1  E-10 

1  .  0E-C9 

3o 

2.2 

1  .8 

5.0 

1  .2 

3-2 

2.4 

2.0 

5.6 

1  .4 

3.1 

2.7 

2.2 

6.2 

1  .5 

3.1 

2.3 

2.3 

6.6 

1  .6 

3.1 

3.0 

2.5 

7.0 

1  .7 

3-1 

3-0 

2.6 

7.5 

1.3 

3.1 

3.1 

2.7 

3.0 

1  .9 

3.1 

3.2 

2.3 

8.3 

2.0 

E+ 1  7 

3.1 

3.2 

2.9 

3.7 

2.0 

3.1 

3-4 

3-4 

1 . 1 E-09 

2.5 

CM  OO  .-t  t(MO  t'- ro  C7>  »  CM  OO  --J  m  VO  t 


Saha  Threshold 

Values  (  /cm  ) 

. )  All  states 

N(i)  (i>7) 

N+(p)  (  p  >7  ) 

2.2E>^  3 

4.5S+1 7 

6 . 4  E  -  T  7 

2 . 5  S  + 1  3 

3 . 6E  +  1  7 

6.3E+17 

2 . 9E+1 3 

3 . 0  E+ 1  7 

7 . 5E+1  7 

3  .  -‘H>1  3 

2 . 4  E+ 1 7 

7 . 8E+1  7 

*  3  •  9  E  +  ’  3 

*  2.0E+17 

3.2E+1 7 

ding  MC3) 

2  .  5 E - J  E  5.CE-G5 

’• w: :CJ 
F igure  1 . 


■  Population  densities  for  M  (sol  LG  L n tE  y  top  3  l. 

^  1 

t=0) ,  N  (solid  line,  bottom  at  t=0) ,  N 
dashed  line),  and  Me  (same  as  M  ).  Ail  figures 
•  are  for  Te=2,0  e.v. 
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CD 


Population  densities  for  N  uv  emitting  states; 
N+(4)  (solid  line,  top),  N+(5)  (dotted  line, 
top),  N  (6)  (dashed  line,  top),  N_r(7)  (solid 

-L 

line,  middle),  N(8)  (dotted  line,  bottom), 

N  (9)  (dashed  line,  bottom),  and  N+(12)  (solid 
line,  bottom) . 


Figure  9, 

Collisional-radiative  recombination  coefficient 
for  optically  thin  neutral  (solid  line)  and  ion 
(dotted  line)  cases. 


.2E-QC  2.5E-0S  S.CE-OS 

TI ME (5) 

Figure  10. 

Collisional-radiative  ionization  coefficients  for 
the  optically  thin  neutral  (solid  line)  and  ion 
('dotted  line)  cases. 
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